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ABSTRACT 

Aims. We present a new sample of X-ray selected galaxy groups and clusters serendipitously observed with the X-ray Telescope 
(XRT) on board of the Swift satellite. Using the XRT archive as of April 2010, we searched for extended sources among 336 GRB 
fields with galactic latitude \b\ >20°. Our selection algorithm provides us with a flux-limited sample of 72 X-ray groups and clusters 
with a well defined selection function and an expected negligible contamination. The sky coverage of the survey goes from the total 
40 deg 2 to 1 deg 2 at a flux limit of 10~ 14 erg s cm -2 (0.5 - 2 keV). This paper provides a description of the XRT data processing, 
the statistical calibration of the survey, and the catalog of detected cluster candidates. 

Methods. All the X-ray sources are detected in the Swift-XRT soft (0.5 - 2 keV) images with the algorithm wavdetect. A size 
parameter defined as the half power radius (HPR) measured inside a box of 45x45 arcsec, is assigned to each source. We select 
extended sources by applying a threshold on the HPR. Thanks to extensive simulations, we are able to calibrate the threshold value, 
(— | _ which depends on the measured net counts inside the box and on the image background, in order to identify all the sources with 

a probability ^ 99% of being extended. The net counts associated to each extended source are then computed by simple aperture 
I ■ photometry. 

Results. We compute the logN-logS of our sample, finding very good agreement with previous deep cluster surveys. We did not 
find any correlation between the cluster and the GRB positions. A cross correlation with published X-ray catalogs shows that only 9 
sources were already detected, none of them as extended. Therefore, ~ 90% of our sources are new X-ray detections. We also cross 
correlated our sources with optical catalogs, finding 20 previously identified clusters. Overall, about ~ 65% of our sources are new 
detections, both as X-ray sources and as clusters of galaxies. 

Conclusions. The XRT follow-up observation of GRBs is providing an excellent serendipitous survey for groups and clusters of 
galaxies, mainly thanks to the low background of XRT and its constant angular resolution across the field of view. A significant 
fraction of the sample (~ 33%) has spectroscopic or photometric redshift thanks to a cross-correlation with public optical surveys. 

Key words, surveys - catalogs - galaxies: clusters: general - galaxies: high-redshift - cosmology: observations - X-ray: galaxies: 
C^l 1 clusters - surveys 
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. 1. Introduction surveys are based on source samples selected by ROSAT, and 

confirmed through optical imaging and spectroscopic observa- 
■ X-ray observations of clusters of galaxies over a significant tions The mQst recent constraints on cosmo i ogica i parameters 
L; ; range of redshifts have been used to investigate the chemical and from x _ ray dusters are bafjed on the chandra f n ow _ up of 400 
. ^ . thermodynamics evolution of the X-ray emitting Intra Cluster deg 2 rqSAT serendipitous survey and of the All-Sky Survey 
^ ■ Medium (ICM, see Etton et al. 2004; Balestra et al. 2007; (V ikhlinin et al. 2009; Mantz et al. 2010). Renewed interest in 
»H . Maughan et al. 2008; Anderson et al. 2009), and to constrain the fieM of cosmological tests with clusters has been recently 
w . the cosmological parameters and the spectrum of the primordial provided by Sunyaev-Zel'dovich (SZ) surveys from the South 
density fluctuations (Rosati et al. 2002a; Schuecker 2005; Voit p ole Telescope Survey (R e ichardt et al. 2012) and the Atacama 
2005; Borgam 2008; Vikhlimn et al. 2009; Mantz et al. 2010; Cosmology Pro j e ct (Sehgal et al. 201 1). At present, only mod- 
Allen et al. 201 1). In this respect, X-ray surveys of clusters of est i m p rovem ents have been obtained with respect to constraints 
galaxies represent a key tool for cosmology and the physics of from WMAP7 plus baryonic acoustic oscillations plus Type la 
large scale structure. The need of assembling larger and larger supernova . The future of sz cluster surveys is very prom i s i ngi 
X-ray selected cluster samples with well denned completeness but at present an x _ ray follow _ up of sz clusters is still neededi 
criteria is one of the critical issues of present-day cosmology. In eimer for narrowing down the uncertainties on the cluster mass, 
order to build statistically complete cluster catalogs, a wide and or tQ firmly evaluate purky and comple teness of the sample. For 
deep coverage of the X-ray sky is mandatory. example, a large effort is being invested in the X-ray follow-up 
To date, there is a remarkable lack of recent wide area X- with Chandra (PI B. Benson) of SZ selected clusters from the 
ray surveys suitable to this scope. Most of the existing cluster South Pole Telescope Survey. This follow-up will provide the 
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X-ray data for 80 massive clusters spread over 2000 deg 2 in the 
redshift range 0.4< z <1.2. 



New X-ray surveys of clusters of galaxies in the Chandra 
and XMM-Newton era are based on the compilation of serendip- 
itous medium and deep-exposure extragalactic pointings not as- 
sociated to previously known X-ray clusters. Among these, one 
of the first survey was assembled by Boschin (2002), who found 
36 clusters (among them 28 new detections) in 5.55 deg 2 of 
surveyed area. Eventually, the Chandra Multiwavelength Project 
(ChaMP) Serendipitous Galaxy Cluster Survey (Barkhouse et al. 
2006) identified about 50 cluster and group candidates from 130 
archival Chandra pointings covering 13 deg 2 . Of the 50 clusters, 
about 16 are expected to have redshift z >0.5. 

More effort is devoted to the XMM-Newton archival data, 
also motivated by the larger solid angle and the nominal larger 
sensitivity. The XMM-Newton Distant Cluster Project (XDCP, 
Fassbender et al. 2011) take advantage of an optical and IR 
follow-up of extended source candidate in XMM images, to 
identify high-z clusters. So far, the XDCP survey has yielded 
22 spectroscopically confirmed clusters in the redshift range 
0.9< z ~ 1.6. A key step in XDCP is the identification of high-z 
candidates based on optical and NIR photometric technique, as 
well as extensive spectroscopic work. In this respect, it cur- 
rently provides the largest sample of confirmed galaxy clus- 
ters at z >0.8, and the purity of the sample is extremely high. 
Another small survey (2.1 deg 2 ) conducted using XMM archival 
data is SEXCLAS (Kolokotronis et al. 2006), which include 19 
serendipitous detections down to 6 x 10~ 15 erg s _1 cirT 2 . 

The largest project based on the entire XMM archive 
is the XMM Clusters Survey (XCS, Romeretal. 2001; 
Lloyd-Davies et al. 2011), with about 1000 cluster candidates 
over a solid angle of 410 deg 2 (the largest based on current 
X-ray facilities). Recently, about 500 clusters have been opti- 
cally identified out of ~ 1000 candidates (Mehrtens et al. 2012). 
Another survey based on the XMM archive is XCLASS, which 
is limited to brighter sources and aims at constraining cosmo- 
logical parameters on the basis of the X-ray information only 
(Clerc et al. 2012). Finally, a survey combining 41.2 deg 2 of 
XMM and Chandra overlapping archival data, plus 122.2 deg 2 of 
Chandra only, has been presented in Peterson et al. (2009), for a 
total of 462 new serendipitous sources. Clearly, there is a signif- 
icant overlap among the cluster samples derived from serendipi- 
tous XMM/Chandra surveys. 

In addition, there are also ongoing dedicated, contiguous 
surveys. The cluster survey in the Subaru-XMM Deep Field 
(SXDF) reaches a depth of 2 x 10 -15 erg s _1 cm -2 over 1.3 
deg 2 with 57 X-ray clusters identified with the red-sequence 
technique (Finoguenov et al. 2010). Similar depth has been 
reached in the COSMOS field over 2 deg 2 , for a total of 72 
clusters (Finoguenov et al. 2007). The XMM-Newton-Blanco 
Cosmology Survey project (XMM-BCS, Suhada et al. 2012) is 
a multiwavelength X-ray, optical and mid-infrared cluster sur- 
vey covered also by the South Pole Telescope and the Atacama 
Cosmology Telescope with the aim of studying the cluster pop- 
ulation in a 14 deg 2 field. The analysis of the first 6 deg 2 pro- 
vided a sample of 46 clusters. Finally, the largest contiguous 
survey is the XMM Large Scale Structure Survey (XMM-LSS, 
Pierre et al. 2007) which is covering a region of 1 1 deg 2 with the 
aim of tracing the large scale structure of the Universe out to a 
redshift of z ~1. At present, a first data release from an area of 
6 deg 2 consists in 66 spectroscopically confirmed clusters with 
0.05<z <1.5 (Adamietal. 2011). 



One of the most important goal of these surveys is to find 
massive clusters at high z. Clearly, tracing the most massive, 
gravitationally bound structures in the Universe up to the highest 
possible redshift is of paramount importance for structure forma- 
tion and for cosmology. While the redshift limit for X-ray se- 
lected clusters is z ~ 1 .57 (in XDCP), extended X-ray emission 
has been detected in optical and IR selected clusters at z = 1 .75 
(Stanford et al. 2012), z = 2.07 (Gobat et al. 201 1), with an ex- 
treme candidate at z ~ 2.2 (see Andreon & Huertas-Company 
2011). 

The present situation is summarized in Table 1. This pic- 
ture is not expected to change significantly in the next years, 
with a modest increase in the source statistics and in the qual- 
ity of the X-ray data. In particular, the study of distant (z > 1) 
clusters relies almost exclusively on time-expensive follow-up 
with Chandra. In the case of Chandra, the process of assembling 
a wide and deep survey is slow due to the small field of view 
(FOV) and the low collecting area. In the case of XMM-Newton 
the collecting area is significantly higher, but the identification 
of extended sources, in particular at medium and high redshift, is 
more difficult due to the larger size of the Point Spread Function 
(PSF, whose half energy width is 15" at the aimpoint) and its 
degradation as a function of the off-axis angle. This is not an 
issue for nearby or medium redshift clusters and groups, but it 
becomes a problem at z > 1, where the optical and IR data play a 
dominat role in cluster identifications. In addition, the relatively 
high and unstable background may hamper the proper character- 
ization of low surface brightness sources. In conclusion, while 
their design is optimized to obtain detailed images of isolated 
sources to explore the deep X-ray sky, a substantially different 
mission strategy is required for surveys. 

In this work we present a new X-ray cluster survey using 
the rich archive of the X-ray Telescope (XRT) onboard of the 
Swift satellite (Burrows et al. 2005) which has never been used 
for the detection of extended sources. Despite its low collect- 
ing area (about one fifth of that of Chandra at 1 keV, see Figure 
1), XRT has two characteristics which are optimal for X-ray 
cluster surveys: a low background and an almost constant PSF 
across the FOV (Moretti et al. 2007). Moreover, almost all the 
Swift/XRT pointings can be used to build a serendipitous survey. 
In fact, the operation mode of the XRT observations considered 
in this work, consists of a prompt follow-up of fields centered on 
Gamma Ray Bursts (GRB) detected by Swift. As we will show 
in § 6, the GRBs do not show any correlation with our extended 
sources. Therefore, despite its small size, the XRT can be suc- 
cessfully used to build an unbiased cluster survey. 

The catalog of extended sources presented in this paper con- 
stitutes the Swift X-ray Cluster Survey (SXCS) and it is based 
on the 336 GRB fields with galaxtic latitude \b\ >20° present 
in the XRT archive as of April 2010. The sources are identified 
thanks to a very simple but effective algorithm to select extended 
X-ray sources in the XRT soft band images. In this paper we 
adopt a conservative detection threshold which guarantees very 
low contamination and a well defined completeness function. 
This catalog is complementary to the catalog of point sources 
identified in the GRB fields of XRT by Puccetti et al. (2011), 
which focus on the study of the AGN population. 

The paper is organized as follows. In § 2 we provide a de- 
scription of the principal characteristics of the XRT. In § 3 we 
describe the field selection and the data reduction. In § 4 we de- 
scribe the detection algorithm and the selection of the extended 
source candidates. In § 5 we describe how we performed ex- 
tended source photometry. In § 6 we present the final list of 
groups and clusters in our sample, compute the sky coverage 
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Table 1 : Flux limited X-ray Cluster Surveys 



Name 


Flux limit 


solid angle 


Number of sources 


Reference 




cgs (0.5 - 2 keV) 


deg 2 






SEXCLAS 


0.6 x 1(T 14 (min) 


2.1 


19 


Kolokotronis et al. (2006) 


DCS 


0.6 x 10- 14 (min) 


5.55 


36 


Boschin (2002) 


ChaMP 


1.0 x 10" 14 (min) 


13.0 


49 


Barkhause et al. (2006) 


SXCS 


1.0 x 10-' 4 (min) 


40.0 


72 


This work 


XDCP 


1.0 x 10~ 14 (average) 


76.0 


22 (z > 0.9) 


Fassbender et al. (201 1) 


XCLASS 


2x 10" 14 (min) 


90.0 


347 


Clerc et al. (2012) 


Peterson09 


~ 0.3 x 10- 14 (min) 


163.4 


462 


Peterson et al. (2009) 


XCS 


> 300 net cts 


410.0 


993 


Lloyd-Davies et al. (201 1) 


SXDF 


0.2 x 10~ 14 (min) 


1.3 


57 


Finoguenov et al. (2010) 


COSMOS 


0.2 x 10~ 14 (min) 


2.1 


72 


Finoguenov et al. (2007) 


XMM-BCS 


0.6 x 10- 14 (min) 


6.0 


46 


Suhada et al. (2012) 


XMM-LSS 


~ 10- 14 (min) 


11.0 


66 


Adamietal. (2011) 



Notes. List of X-ray flux limited cluster surveys based on Chandra or XMM data, updated at the time of writing (May 2012), plus the Swift/XRT 
Cluster Survey presented in this work (highlighted in bold). Surveys based on archival data are listed in the upper part of the table, while dedicated 
(contiguous) surveys are shown in the lower part. Surveys are ranked according to the total solid angle. The total number of clusters refer to the 
X-ray selected only, while the quoted solid angle is the maximum covered by the survey. Note that the limiting soft fluxes quoted in this table 
(when available) may refer to the lowest value in the sample (minimum) or to an average value over the entire solid angle (average). For a full 
characterization of the survey depth, i.e., the sky coverage as a function of the flux, we refer to the corresponding papers. 




0.1 1.0 10.0 

E (keV) 



Fig. 1: Effective area of XRT (green solid line) as a function of 
the energy compared with that of the Chandra satellite (red dot- 
ted and cyan dashed lines for Chandra ACIS-S and ACIS-I, re- 
spectively). 



and the LogN-LogS, and check for possible selection bias. In 
§ 7 we correlate our catalog with existing databases to identify 
previously known sources and collect the spectroscopic or pho- 
tometric redshifts of member galaxies available in the literature. 
In § 8 we discuss our results in the context of current and future 
X-ray surveys. Finally, in § 9 we summarize our conclusions. 



2. XRT characteristics 

The XRT is part of the scientific payload of the Swift satel- 
lite (Gehrels et al. 2004), a mission dedicated to the study of 
gamma-ray bursts (GRBs) and their afterglows operating since 
January 2005 '. GRBs are detected and localized by the Burst 
Alert Telescope (BAT, Barthelmy et al. 2005), in the 15-300 
keV energy band and followed-up at X-ray energies (0.3-10 
keV) by the X-ray Telescope. The XRT (Burrows et al. 2005) 
is an X-ray CCD imaging spectrometer which utilizes the third 
flight mirror module originally developed for the JET-X program 
(Citterio et al. 1994). The mirror module focuses X-rays (0.2-10 
keV) onto a XMM-Newton/EPICMOS CCD detector consisting 
of 600 x 600 pixels, with a nominal plate scale of 2.36 arcsec - 
onds per pixel, which provides an effective FOV of the system 
of ~ 24 arcmin. The PSF, similar to XMM-Newton, is character- 
ized by an half energy width (HEW) of ~ 1 8 arcsec at 1 .5 ke V 
Moretti et al. (2007). The PSF dependence on the off-axis angle 
is very weak. This is due to the fact that the CCD is intentionally 
slightly offset along the optical axis from the best on-axis focus 
in order to have a uniform PSF over a large fraction of the FOV, 
with a negligible dependence on the photon energy. To show the 
remarkably flat behaviour of the PSF, we show in Figure 2 the 
measured Half Power Radius (HPR) within a box of 45 x 45 arc- 
sec as a function of the off-axis angle for all the sources detected 
in the XRT fields used in this work. The PSF can be analytically 
described by a King function with slope /3 ~ 1 .45 and core radius 
r c ~ 5.3 arcsec (Moretti et al. 2007). 

One of the most relevant aspect for our purposes is the low 
level and the high reproducibility of the background not asso- 
ciated to astronomical sources (NXB). Due to the low orbit and 
short focal length, the NXB is the lowest among the currently op- 
erating X-ray telescopes (see Hall et al. 2007). As a reference, 
the NXB of XRT calculated per solid angle and normalized by 



1 In 2012 the NASA Senior Review committee allocated full funding 
for for Swift in the period 2013-14 and recommended full funding also 
for 2015-16 with next review in 2014. 
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Fig. 2: Measured HPR within a box of 45 x 45 arcsec as a func- 
tion of the off-axis angle for all the sources detected in the XRT 
fields used in this work. 



500 Swift Gamma-Ray Bursts (Jan05 -Apr10) 




Fig. 3: Distribution on the sky (Aitoff projection) of the GRB 
fields in the XRT archive as of April 2010 (red dots). The 
fields selected in our survey are those outside the Galactic plane 
(shaded area). 

the effective area is a factor ~ 7 lower than Chandra in the 0.5- 
7.0 keV energy band (Moretti et al. 2012 in preparation). 

3. Fields selection and data reduction 

We consider the entire Swift/XRT archive from February 2005 
to April 2010, including 502 fields in the corresponding GRB 
positions, distributed in the sky as shown in Figure 3. We select 
the fields with Galactic latitude \b\ >20°, to avoid crowded fields 
and strong Galactic absorption. This selection provides us with 
336 fields which we use to search for extended sources. 

In this paper we use photon counting mode data with the 
standard grade selection (grade 0-12) reduced by means of the 
xrtpipeline task of the current release of the HEADAS soft- 
ware (version v6.8) with the most updated calibration at the time 
of writing (CALDB version 20 1 1 1 03 1 , Nov 20 1 1 ) . More detail 
on the standard data reduction can be found in the instrument 



user guide 2 . Then, we proceed with a customized data reduc- 
tion, aimed at optimizing our data to the detection of extended 
sources. First we exclude the external CCD columns (corre- 
sponding to the detector coordinates Detx < 90 and Detx > 510) 
which are affected by the presence of out-of-time-events from 
corner calibration sources (see Moretti et al. 2009, for a detailed 
map of the CCD and a discussion on the XRT background). The 
corrected FOV, after removal of the external CCD columns, is 
16.5' x 18.9' (0.087 deg 2 ) for each pointing. Different obser- 
vations of the same object and corresponding exposure maps are 
merged by means of the the extractor and farith tasks of the 
HEADAS software, respectively. 

In order to further reduce the background, we investigate the 
background light curves in the soft image of each field. At the 
beginning of each observation, corresponding to the maximum 
emission from the GRB, the background on the entire image is 
significantly larger than the typical, quiescent value. Actually it 
was previosly known (Moretti et al. 2009) that for very high flu- 
ences 3 , about 5% of the GRB X-ray emission is scattered across 
the image. 

We effectively tested that by removing all the data taken be- 
fore the epoch To + 0. 1 x t exp , where 7*o is the epoch of the begin- 
ning of the observation, and t exp is the total effective exposure 
time, the average background on a typical image can be reduced 
by ~ 5%. This empiric rule is motivated by the fact that the to- 
tal exposure time of a typical GRB follow-up observation is set 
by the GRB emission itself, therefore brighter GRBs will have 
larger observation length and therefore larger time intervals re- 
moved. In addition, note that the time interval 0.1 x t exp does 
not correspond to 10% of the exposure time, since this is dis- 
tributed over several orbits, each one with an observability time 
of 1500 s over a total orbit period of 5400 s. Figure 4 shows an 
example of a typical light curve, where the time interval that has 
been removed is shown in grey. The removal of the data up to 
the time 7o + 0.1 x t exp corresponds typically to the first two or 
three orbits. On average, the time interval 0.1 x t exp corresponds 
to 2-3% of the effective exposure time, since this is distributed 
over several days with 10-20 ks observed every day. We also note 
that removing the first orbits has the effect of reducing the non- 
homogeneity of the background caused by the GRB emission, 
which would have affected the detection of extended sources. In 
Figure 5 we show the distribution of the effective exposure times 
of the XRT fields after this correction. 



4. Selection of extended-source candidates 

4. 1 . Initial source detection 

The identification of the X-ray sources is performed with 
the standard algorithm wavdetect within the CI AO soft- 
ware used successfully on Chandra and XMM-Newton images 
(Lloyd-Davies et al. 2011). We run the algorithm on the im- 
ages obtained in the soft (0.5 - 2 keV) band using a set of 
six wavelet scales corresponding to 2, 4, 8, 16 and 32 image 
pixels. We use a mild selection threshold corresponding to the 
wavdetect threshold parameter 10~ 6 . Despite not guarantee- 
ing a high completeness nor high purity, this choice is moti- 
vated by the fact that the selection of extended sources will be 
based on the growth curve method, and that only relatively high 
S/N sources will be selected as we will show in Section 4.5. 



2 http://heasarc.nasa.gov/docs/swift/analysis/documentation 

3 The fluence is the total energy delivered per unit area, obtained by 
integrating the source flux over time. 
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Fig. 4: Light curve of the background of the field of GRB051008 
plotted as a function of the time elapsed from the GRB trigger. 
The saw-tooth appearance is due to the orbits of Swift (one orbit 
lasts 90 minutes). The removed time interval below the epoch 
To + 0.1 x t exp is shown in grey. 
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Fig. 5: Distribution of the exposure times in ks for the 336 se- 
lected GRB fields. 



Therefore the wavdetect threshold parameter has no effects on 
our final cluster candidates list. 

The soft band is optimal to identify extended ICM emission, 
since in this band XRT has the highest effective area and the 
lowest instrumental background. In this way we maximise our 
sensitivity in the detection of extended source powered by ther- 
mal bremsstrahlung, whose soft band emission has a small K- 
correction up to z ~ 1, at least for hot clusters (kT>3 keV). 
We also considered modifications to the standard soft band in 
order to find the energy range optimal for cluster detection, sim- 
ilarly to what has been done for Chandra and XMM-Newton 
(Scharf 2002). We explored the use of a narrower energy band 
in order to reduce the effect of the Galactic background, which 
rapidly grows below 1 keV. We found that a narrower energy 
band would have a minor positive effect on the detection of clus- 
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£ 0.0 



-0.5 
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Fig. 6: Hardness ratios modeled for XRT as a function of redshift 
for groups and clusters with several temperatures (red continu- 
ous lines, kT respectively equal to 0.5, 1.0, 2.0, 3.0, 5.0, 7.0, 10.0 
and 12.0 keV from the bottom to the top) and for AGN with dif- 
ferent intrinsic absorption (blue dashed lines, N# respectively 
equal to 0.01, 0.1, 0.3, 1.0, 3.0 and lO.OxlO 22 cm" 2 from the 
bottom to the top). In the top right, we show the typical error bar 
on HR for a source with 200 total net counts. 



ters, while negatively affecting the detection of low temperature 
groups whose emission is below 2 keV. We also explored the use 
of a larger upper energy value, but the poor photons statistics for 
thermal spectra at high energies and the rapidly increasing hard 
background nullify any advantage for values above 2 keV. Since 
the small, positive effects we found in changing the energy band 
are also significantly dependent on the temperature of the ICM, 
we decide to maintain the standard 0.5 - 2 keV band. 

Incidentally, we notice that, given the sensitivity at high en- 
ergies of the XRT, it is possible in principle to use also the hard 
band (2-7 keV) to disentangle thermal emission from power-law 
emission typical of Active Galactic Nuclei (AGN), on the ba- 
sis of the source hardness ratio. The hardness ratio is defined 
as HR=(H-S)/(H+S), where H and S are the hard band and soft 
band counts respectively, corrected for vignetting. In Figure 6 
we show the hardness ratio as a function of redshift expected for 
groups and clusters at different temperatures, compared to the 
hardness ratio of AGN with different intrinsic absorption N#. 
Values for clusters and AGN are taken using Xspec vl2.6.0 
assuming in the first case an absorbed thermal model with a 
mean Galactic absorption of Nh = 3 x 10 20 crrT 2 , and in the 
second case assuming an absorbed power law with an intrin- 
sic redshifted absorption. We find that only nearby, strongly ab- 
sorbed AGN have hardness ratio values clearly different from 
those of the ICM. We conclude that the hard band does not pro- 
vide a relevant information for the detection of groups and clus- 
ters. Therefore we will use only the soft band in this work. The 
full band will be used for the spectral analysis of the brightest 
sources in a companion paper (Moretti et al. in preparation). 
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Since the algorithm wavdetect has never been tested thor- 
oughly on XRT images, we cannot immediately associate an 
expected contamination level to the threshold value of 10 6 . 
In any case, at this stage we are mostly interested in collect- 
ing all the possible extended-source candidates, while the con- 
tamination level will be estimated during our selection process. 
This step leads us to identify a total of ~ 10 4 sources in the 
336 GRB fields. This number is somewhat lower than that ex- 
pected on the basis of the Swift serendipitous survey (SwiftFT, 
see Puccetti et al. 2011) which was based on a ~ 20% smaller 
field selection. The lowest number of net photons in our total 
source list is as low as 5. Despite the low background, Poisson 
fluctuations may appear as spurious sources at such a low sig- 
nal. However, we do not apply any filter to this parent list, since 
much more stringent thresholds will be applied when selecting 
extended sources, as described in the next subsection. 

4.2. Source extent determination 

For each source we consider a 60x60 pixels box centered around 
the source. In each of these boxes, we mask the other sources 
whose emission can overlap with the central one, by removing 
a PSF image. The PSF is accurately fitted with a King profile 
normalized to each source photometry assuming a core radius of 
r c = 5.3 arcsec and a fi = 1.45. These values are valid at en- 
ergies ~ 1.5 keV and have been recomputed in-flight by fitting 
the many high S/N sources observed by XRT, therefore updat- 
ing the values in Moretti et al. (2005) and Moretti et al. (2007). 
This gives us a sub-image cleaned from point-like sources, with 
the considered source at its center. In order to account for defects 
due to missing columns or removed pixels in the CCD, we divide 
this image by the corresponding soft exposure map. Then, we 
compute the growth curve of the source within a box of 45 x45 
arcsec. This region corresponds to 2.5 times the HEW (which is 
18 arcsec) and includes 80% of the flux of a point source. From 
the growth curve within this region, we measure the HPR, effec- 
tively defined as the 50% encircled energy radius within a box 
of 45 x 45 arcsec. Note that the HPR is smaller than the HEW, 
since it refers to only 80% of the total flux. The choice of a re- 
stricted regions is motivated by the need of sampling the growth 
curve with a good signal-to-noise ratio (S/N) for the majority of 
the sources in the initial list. This choice does not affect by any 
means the final total flux, which will be computed according to 
a different procedure as shown in the next subsections, and it is 
used only to classify extended versus unresolved sources. 

4.3. Point source simulation, and determination of threshold 
parameters 

Extensive simulations spanning all the parameter range found 
in the survey (in particular source fluxes and background lev- 
els) have been used to estimate the expected distribution of the 
HPR for point sources. The input flux distribution for unre- 
solved sources is consistent with the number counts obtained 
by Puccetti et al. (2011) and with the deeper CDFS counts 
(Rosati et al. 2002b), while the fluxes for clusters are distributed 
according with the number counts obtained in the ROSAT Deep 
Cluster Survey Rosati et al. (2002a). 

The distribution of the measured HPR in the simulations as 
a function of the source counts within the 45x45 arcsec box 
is shown in Figure 7. The simulations (including only point 
sources) are used to define the threshold values HPR,;, above 
which a source is inconsistent with being unresolved at the 99% 
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Fig. 7: Distribution of the measured HPR within a 45x45 arc- 
sec box for simulated unresolved sources as a function of the 
soft counts within the box. The same range of exposure times 
and backgrounds present in the real data has been adopted in 
the simulations. The continuous line shows the 99% confidence 
level envelope of the source distribution. 



level (see red line in Figure 7). This criterion is extremely sim- 
ple and it does not depend on the off-axis angle 0, thanks to 
the flat PSF. However, HPR,;, is a function of the source counts 
measured within the 45x45 arcsec box, in particular it increases 
significantly below 60 counts. Above this value, HPR,;, ranges 
between 3 and 4 image pixels (equal to 7.1-9.4 arcsec for a pixel 
size of 2.36 arcsec). Instead of using a single HPR,;, based on 
the entire simulation, we compute HPR,;, for a set of different 
background values, finding a significant dependence which will 
be taken into account in the source selection. By directly ap- 
plying this criterion, we expect to have a number of spurious 
sources equal to 1 % of the total source number. We remark that 
the counts within the 45x45 arcsec box can be much lower than 
the total net counts, particularly for extended sources. 

This procedure is applied to all sources detected with 
wavdetect. Figure 8 shows the distribution of the HPR of 
the detected point sources (grey dots) versus the source counts 
measured inside the HPR, with the extended source candidates 
marked as black and red big dots. Note that extended source can- 
didates are selected on the basis of a HPR,;, value which depends 
also on the background (tied to the exposure of the image). The 
HPR,;, for different exposures are shown in Figure 8 as solid lines 
ranging from 1 ks (orange line) to 1 Ms (blue line). The number 
of extended source candidates grows rapidly when the number 
of detected counts within the box decreases towards low values. 
The inclusion of all the candidates in our survey clearly would 
imply a large number of spurious sources, given the large total 
number of sources (~ 10 4 ). We finally apply the threshold on 
the total number of counts to be > 100 in the extraction radius 
(defined as the circular region where the source surface bright- 
ness is larger than the background level). The red big dots indi- 
cate the candidate extended sources matching all our selection 
thresholds. 

The sharp threshold on the soft net counts allows us to com- 
pute directly a sky coverage as a function of the source flux 
within the extraction region (see Section 6). Clearly, the sky cov- 
erage information would be sufficient to compute the logN-logS 
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only in absence of any morphology bias. Actually, we do miss 
a fraction of extended sources, mostly due to the wide range of 
intrinsic morphologies. The completeness of our survey is then 
estimated as described in the next subsection. 

4.4. Completeness and contamination level 

To properly characterize the quality of our sample we need to 
assess the completeness, defined as the fraction of extended 
sources actually selected as extended by our procedure as a func- 
tion of their soft net counts. In order to estimate the complete- 
ness, we realize another set of simulations with a different strat- 
egy. We simulated a few thousand extended sources, spanning 
a wide range of fluxes consistent with the number counts mea- 
sured in the ROSAT Deep Cluster Survey (Rosati et al. 1998). 

Our extended sources are modelled starting from ten real 
cluster images originally obtained with the Chandra satellite 
(therefore with a resolution much higher that that of XRT), 
cloned at a typical redshift and resampled at the XRT reso- 
lution. This cloning procedure, already used to investigate the 
evolution of cool core clusters at high redshift (see Santos et al. 
2008, 2010) allows us to measure our completeness for a real- 
istic range of surface brightness profiles, representing the local 
population of groups and clusters of galaxies: temperatures go 
from 2 keV to 8 keV, the redshift of the clusters are all below 0.2, 
and half of the clusters have a cool core. Clearly, the limited set 
of group/cluster templates, and the assumption of no evolution 
in the surface brightness properties with redshift, may introduce 
some difference with respect to the actual cluster population. The 
impact of the intrinsic morphologies on the source selection is 
properly taken into account as long as the morphologies of the 
simulated sources are representative of the real sources. This is 
admittedly a limitation of the present approach. Still, this proce- 
dure is the most accurate with the present knowledge. The sim- 
ulated clusters are positioned randomly on real XRT images and 
the fluxes are converted into count rate using the response func- 
tion of the instrument, including vignetting effects. We run our 
detection procedure on the simulated images, and compute the 
completeness as the recovered source fraction (i.e., the fraction 
of sources detected and characterized as extended with our crite- 
rion) as a function of the input counts. As shown in Figure 9, we 
reach a completeness level ~70% for sources above our thresh- 
old of ~100 net counts within the extraction regions. This level 
increases to ~9Q% for sources with ~200 net counts, and reach 
a completeness level of >95% for sources with ~300 net counts. 
We will use this information to correct the logN-logS of our sam- 
ple in § 6. We remark that this completeness function is by no 
means general: it depends, in fact, on the intrinsic properties of 
the survey, on the exposure time distribution, and on the source 
selection and classification method. We are actually aiming at in- 
creasing the completeness down to a threshold of ~50 net counts, 
while maintaining a low contamination level, thanks to a differ- 
ent detection algorithm based on Voronoi tessellation (Liu et al. 
in preparation). 

Finally, our simulations allow us to estimate also the con- 
tamination level. This is obtained directly by averaging the num- 
ber of spurious sources surviving our selection criteria in all our 
survey realizations. The average expected number of spurious 
sources in the SXCS survey turns out to be ~ 5. The contamina- 
tion is due to Poisson noise in the measured parameters (HPR, 
photometry) and to blending. We neglect any spatial correlation 
which, however, is expected to give a negligible enhancement to 
the contamination level. 
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Fig. 8: Distribution of the measured HPR within a box of 45 x45 
arcsec for all the sources detected in the XRT fields as a function 
of the soft counts measured within the HPR. All possible ex- 
tended source candidates are marked with black big dots, while 
red big dots indicate sources with >100 soft net counts from the 
aperture photometry in the extraction radius (see Section 5) and 
are therefore included in the cluster candidate list. The contin- 
uous lines show the typical 99% confidence level envelope for 
different exposures, from 1 ks (orange line) to 1 Ms (blue line). 



4.5. Final sample selection 

Our algorithm identified 87 candidate extended sources with 
>100 soft net counts. We proceeded to a careful visual inspec- 
tion of these candidates, including a look through the optical im- 
ages taken from the DSS and from the SDSS when available. 
This visual inspection allows us to remove 3 sources which are 
identified with stars 4 . In addition, 3 sources are removed since 
they are identified with nearby large galaxies. Note that all these 
sources are correctly identified as extended by our algorithm. 

We also remove 6 sources since they are easily recognized 
as blended by visual inspection. This effect was included in the 
simulations which predicted a number of spurious sources equal 
to 5. We argue that our visual inspection actually remove most 
of the spurious detections. The visual inspection is made pos- 
sible thanks to the limited size of our sample. For larger sam- 
ples, a more efficient, self-consistent detection algorithm would 
be needed in order to keep the contamination level under con- 
trol down to lower fluxes, without recurring to time-consuming 
visual inspections. 

Another 3 sources are found to be very close to the edge of 
the XRT images, where the gradient of the exposure map signif- 
icantly affects the photometry. We remove them since their pho- 
tometry would necessarily rely on a substantial extrapolation of 
their properties. The final sample includes 72 extended sources. 
Figure 10 shows a typical example of a GRB field imaged by 
XRT. In this case (GRB050505) three extended sources, high- 
lighted by green circles, have been selected by our algorithm. 



4 Bright stars appear as extended in the XRT image due to the satura- 
tion of the brightest pixels. Despite this, it is easy to remove these spuri- 
ous extended sources after a crosscorrelation with any optical database. 



7 



E. Tundo et al.: The Swift X-ray Telescope Cluster Survey 





1 00 1 000 

soft net counts 



Fig. 9: Completeness level, expressed as the ratio of recovered 
extended sources on the number of simulated extended sources 
as a function of the input soft counts. 



5. Extended sources photometry 

The measured net counts from each source are computed by 
aperture photometry up to an extraction radius R ext , after remov- 
ing the emission associated to point sources included in this re- 
gion. R ex , is defined as the radius where the surface brightness 
profile of the source reaches the background level. The total net 
counts are finally obtained by accounting for the missed flux be- 
yond R ext . In order to measure the lost signal, we fit the sur- 
face brightness SB of every extended source with a King profile 
modeled as SBcc (1 + (r/r c ) 2 ) 3 ^^ 2 (Cavaliere & Fusco-Femiano 
1978) leaving both r c and B free to vary. Then, we extend the 
surface brightness profile up to 2 x R ext . The typical ratio of 
the total estimated flux to that actually measured within R ext is 
Cf = 1.05 - 1.10. Extrapolating the profile to larger radii has a 
modest effect on the final results. We also checked that the errors 
on the r c and B parameters affect the correction factor c/ only at 
a 10% level, and therefore constitute a negligible uncertainty for 
the final source flux. 

The net count rate measured for each source is obtained di- 
viding the net counts by the exposure time, after correcting for 
vignetting effects in the soft (0.5 - 2 keV) band. The average 
vignetting correction in the extraction region is estimated by the 
ratio of the photon-weighted average value of the soft exposure 
map within R e x , to the value of the soft exposure map at the aim- 
point. From the corrected net count rate, the energy flux is com- 
puted simply by multiplying it by the energy conversion factor 
(ECF) computed at the aimpoint for a suitable spectral model. 

The ECF for thermal X-ray emission from the ICM is com- 
puted assuming a mekal model within XSPEC. Since we do 
not know a priori the temperature of our sources, we explore 
a range of temperatures from 1 to 12 keV, and a redshift range 
from to 1 .5. The corresponding ECF values in the soft band are 
shown in Figure 1 1 for a typical Fe abundance of 0.3 Z in units 
of Anders & Grevesse (1989). We note that the ECFs vary less 



Fig. 10: XRT smoothed image of the GRB050505 field in the 
soft (0.5 - 2.0 keV) band; the exposure time is ~170 ks, and the 
image size is ~24x24 arcmin. The three clusters in the fields are 
highlighted by the green cicles. The bright source at the center 
istheGRB. 



than 2% for kT>3 keV at any redshift. The largest changes are 
found for lower temperatures, up to a maximum of 8% for kT-1 
keV. However, given the average L - T relation measured lo- 
cally, and considering that it is approximately constant with red- 
shift (Reichert et al. 201 1; Branchesi et al. 2007; Maughan et al. 
2006; Ettori et al. 2004; Mushotzky & Scharf 1997), we can de- 
rive a minimum temperature detectable at a given redshift for our 
flux limit. If we set Fn m ~ 10~ 14 erg s _1 cirT 2 , we find that clus- 
ters and groups with kT<2 keV can be observed only for z <0.5. 
This allows us to conclude that the maximum variation is be- 
low 4%, and that this number is little affected by the actual ICM 
metallicity. Therefore we assume an average ECF of 2.35xl0~ n 
erg s _1 cirT 2 / (cts s~') with a maximum systematic uncertainty 
of 0.1 x 10~ n erg s _1 cirT 2 / (cts s~'). 

The ECFs, however, have a more significant dependence 
on the Galactic absorption. This can be estimated thanks to 
the Galactic Ng values measured in the radio survey in the 
Leiden/ Argentine/Bonn survey (Kalberla et al. 2005). As shown 
by the histogram in Figure 12, the distribution of Galactic Ng for 
our sources implies changes up to the order of 30%. Therefore 
we use for each detected source the ECF appropriate to the 
corresponding field. We note that self shielding of molecular 
Hydrogen from ambient UV radiation is expected to occur for 
Ng > 5 x 10 20 cm" 2 (Arabadjis & Bre gman 1999), and this 
molecular gas absorbs also X-rays. Therefore for this Ng val- 
ues X-ray fluxes are slightly underestimated. We do not correct 
for this effect (see also Lloyd-Davies et al. 201 1). 

To summarize, the total unabsorbed, soft fluxes, including 
the correction up to 2 R ext , are measured as: 

S .5-2keV = ECF(N H ) X Cf X C rate (1) 
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Fig. 11: Energy Conversion Factor (ECF) as a function of the 
ICM temperature for different redshifts. Here we assumed no 
Galactic absorption. 
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Fig. 12: Energy Conversion Factor (ECF) as a function of the 
Galactic Nh for a typical cluster (kT=5 keV, z =0.4). The his- 
togram shows the distribution of Galactic N# in our cluster sam- 
ple. 



where ECF(Nh) is the energy conversion factor which accounts 
for the Galactic absorption (see Figure 12), cj is the correction 
factor for the flux between 2 x R ext and R ext , and C rate is the soft 
band count rate mesured within R ext and corrected for vignetting 
effects. 



6. Cluster catalog and Number Counts 

Aperture photometry within R ext and total energy flux are shown 
in Table 2 for all the sources of our sample. The net counts er- 
ror is obtained from the Poissonian error on the numbers of net 
photons. We remark that, as explained in detail in the previous 
section, the photometry of our extended sources refers to the ex- 
traction radius R exl defined as the radius where the fitted surface 
brightness falls below the measured local background, while the 
total energy flux include the correction factor c j for the flux lost 
beyond R„,. 

In order to compute the number counts for our cluster sam- 
ple, we need to derive the sky coverage of our survey. The sky 
coverage is determined by the sharp limit in the total net pho- 
tons. For each field we compute a flux-limit map obtained as 
ECF(Nh) x 100/ Expmap{t), where Expmap(t) is the expo- 
sure map in units of effective time, therefore including the ef- 
fect of vignetting. Then, the solid angle covered by the survey 
above a given flux is obtained by measuring the total solid angle 
where the flux-limit is lower than a given flux. In Figure 13 we 
show the resulting sky coverage Q(S) and compare it with the 
sky coverage of the 400 Square Degree ROSAT PSPC Galaxy 
Cluster Survey by Burenin et al. (2007) and with the sky cover- 
age of the ROSAT Deep Cluster Survey (RDCS) by Rosati et al. 
(1998). Unfortunately, sky coverage curves for other on-going 
XMM/Oiandra surveys have not been published. The cumula- 
tive number counts for sources brighter than a given flux S is 
therefore given by 

N(>S) = L CfXSi>s C; 1 /£l(S i ) (2) 

where S is the total flux, S , is the soft flux within R ex , of the i th 
source, c/ is the correction factor for the flux outside R ext and 
Cl(Si) is the solid angle corresponding to 5,-. Finally, C, is the 
completeness factor plotted in Figure 9 which depends on the 
net detected photons of the i' h source. In this way, each source 
is weighted with a factor inversely proportional to the survey 
completeness which depends only on source photometry. 

An alternative procedure to correct for completeness is ob- 
tained by randomly extracting the missing sources and comput- 
ing their average effect. The procedure consists in dividing the 
sources in 3 bins of net photons (100-150, 150-200, 200-300). 
Then, in each bin we randomly add a number of sources drawn 
from a Poissonian distribution centered on the expected number 
of missed sources according to Figure 9. Finally we assign a ran- 
dom exposure and energy conversion factor among those in the 
survey to each mock source, and compute its energy flux. Finally 
the logN-logS is recomputed for each Monte Carlo realization. 
The average logN-logS obtained with this procedure with 10 4 
realizations is in very good agreement with that obtained with 
equation 2, therefore confirming that our treatment of the com- 
pleteness is robust. 

Our logN-logS is shown in Figure 14 with its 68% confi- 
dence limits (shaded red region). The confidence limits are com- 
puted with Monte Carlo realizations. We re-extract the flux of 
each source times assuming its Poissonian uncertainties on the 
net detected counts, including the systematic uncertanty on the 
ECFs, obtaining 10 4 realizations of the logN-logS. Then, at each 
flux, we compute the 68% confidence interval around the mean 
value. 

Our results are consistent with the logN-logS measured with 
the ROSAT Deep Cluster Survey (RDCS, Rosati et al. 1998), 
shown in Figure 14 as the cyan region, down to fluxes ~ 10~ 14 
erg s _1 cirT 2 . We are also consistent with the the logN-logS 
measured in the 400 deg 2 survey (Burenin et al. 2007) and in 
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Fig. 13: Solid angle Q covered by our survey as a function of the 
soft flux. As a comparison, we show also the sky coverage of 
the 400 Square Degree Survey (Burenin et al. 2007, blue dashed 
line) and the sky coverage of the RDCS (Rosati et al. 1998, green 
dotted line). 



Fig. 14: Cumulative number counts, corrected for completeness, 
with 1 <x confidence level for SXCS (red area), compared to the 
logN-logS derived from the ROSAT Deep Cluster Survey (cyan 
area), Rosati et al. (1998). Dot-dashed and dashed lines show 
the fit to the 400d (Burenin et al. 2007) and to the COSMOS 
(Finoguenov et al. 2007) number counts, respectively. 



the COSMOS field (Finoguenov et al. 2007), shown as the dash- 
dotted and dotted lines respectively. 

The good agreement with previous results confirm the proper 
characterization of our sample. Nevertheless, we also perform a 
test against a possible correlation between the GRB position and 
the position of our sources, since the detection of an enhanced 
density of clusters and groups toward the position of GRB would 
bias our survey. In Figure 15 we plot the positions of the sources 
in our sample with respect to the GRB positions (set to 0,0). We 
do not find any hint of an increasing surface density of extended 
sources towards the GRB position. We also tested for such a pos- 
sible bias by measuring the logN-logS for the sources whose 
distance is, respectively, below 400 arcseconds from the GRB 
and between 400 and 800 arcseconds from the GRB. The two 
logN-logS agree with each other within the uncertainties, show- 
ing that there is no correlation between GRB and clusters (see 
also Berger et al. 2007). We can safely assume that our field se- 
lection is unbiased with respect to X-ray clusters. 

In order to check whether we properly treated the effect of 
the Galactic absorption, we split the survey into two segments 
according to the Galactic absorption, above and below Nh = 
3 x 10 20 cnT 2 . We find that the logN-logS computed in the two 
cases are in very good agreement, and therefore we conclude 
that our flux measurements do not appear to suffer any bias from 
Galactic absorption. 

To summarize, the sources included in the SXCS catalog rep- 
resent a sample of group and cluster candidates with a negligible 
contamination, a well defined selection function and a robustly 
estimated completeness, spanning two orders of magnitude in 
flux, and reaching the flux limit of ~ 10~ 14 erg s cmT 2 . 
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Fig. 15: Positions of our sources with respect to the GRB posi- 
tion (set to 0,0). 
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7. Cross-correlation with X-ray and optical 
catalogs. 

We checked for counterparts both in previous X-ray surveys and 
in optical cluster surveys, assuming a search radius of 1 arcmin 
from our X-ray centroid. The results are shown in Table 3. We 
find 9 X-ray counterparts to our sources in ROSAT, ASCA and 
Chandra catalogs, none of them characterized as extended and 
none of them with a published redshift. We also find a total of 20 
previously known, optically identified clusters. Among them, 6 
are found in the Wen+Han+Liu cluster sample (WHL, Wen et al. 
2009), which is an optical catalog of galaxy clusters obtained 
from an adaptive-matched filter finder applied to Sloan Digital 
Sky Survey DR6. We find 4 clusters in the Szabo et al. (2011) 
catalog (AMF) based on a similar method in the SDSS DR6. 

We also report one cluster from the Gaussian Mixture 
Brightest Cluster Galaxy (GMBCG, Hao et al. 2010) based on 
SDSS DR7, which is an extension of the maxBCG cluster cata- 
log (Koester et al. 2007) to redshift beyond 0.3 and on a slightly 
larger sky area (~ 8000 deg 2 ), and one cluster in the MaxBCG 
catalog itself. 

Finally, we find 3 clusters in the Abell catalog (Abell et al. 
1989), two cluster in the Northern Sky optical Cluster Survey 
(NSCS, Galetal. 2003; Lopes et al. 2004), and one clus- 
ter in each of the following catalogs: the Zwicky Cluster 
Catalog (CGCG, Zwicky et al. 1963), the Sloan Digital Sky 
Survey C4 Cluster Catalog (SDSS-C4-DR3, based on DR3, 
Miller et al. 2005; von der Linden et al. 2007), the Edinburgh- 
Durham Southern Galaxy Catalogue (EDCC, Lumsden et al. 
1992). Three out of twenty clusters also have an X-ray coun- 
terpart. 

To increase the number of available redshifts, we also search 
for galaxies with published redshift not associated to previously 
known clusters, within a search radius of 7 arcsec from the X- 
ray centroid of our sources. In Table 3 we report 1 1 galaxies with 
redshift, within a search radius of 7 arcsec, as a complement to 
the redshift of the cluster counterpart. The redshift of the cen- 
tral galaxy candidate is always consistent with the photometric 
redshift of the optical cluster counterpart when present. In the 
4 cases where no optical cluster counterpart is found, we ten- 
tatively assign the galaxy redshift to our X-ray source. We note 
that both cluster and central galaxy counterparts have been found 
by adopting a simple matching criterion. A more refined analy- 
sis of the optical data covering our survey is under way (Tundo 
et al. in preparation). 

To summarize, we have 24 optical redshifts (spectroscopic 
or photometric) published in the literature and associated to our 
cluster candidates. Overall, 46 sources in our catalog are new 
detections, both as X-ray sources and as clusters of galaxies. 
Finally, we remark that a total of 32 SXCS sources fall in SDSS 
fields. Thanks to the SDSS depth, we expect to increase the iden- 
tification of our cluster candidates. In particular we expect to be 
able to identify at least the brightest cluster galaxy up to z < 1 
for those cluster candidates which are not already included in the 
SDSS catalogs. This will also provide an estimate of the photo- 
metric redshift of the host cluster. Figure 16 shows a selection 
of SDSS r-band images of SXCS fields, with X-ray contours 
overlaid in green. 

We stress that we can also attempt a measure of the source 
redshift through the X-ray spectral analysis for a significant 
fraction of the sample. The requirements for a successfull iden- 
tification of the redshifted K a Fe line, shown in Yu et al. (201 1) 
for Chandra, do not apply to most of the SXCS sources. However 
the lower background and the flatter effective area of XRT may 



allow X-ray redshift measurements in a lower S/N regime. A 
preliminary result indicates that the X-ray redshift zx of about 
30% of the SXCS sources can be successfully measured. This 
is an important aspect since the measure of zx will complement 
the photometric and spectroscopic redshifts obtained with opti- 
cal follow-up with the aim of using the cluster sample for cos- 
mological tests. We will present the spectral analysis in a com- 
panion paper (Moretti et al. in preparation). 

8. Future prospects 

SXCS data have different properties from Chandra and XMM- 
Newton data as described in § 2 and § 3, and therefore SXCS 
extended sources have a different selection function. For these 
reasons the SXCS is a valuable touchstone to the cluster cata- 
logs obtained from XMM-Newton and Chandra. The full SXCS 
survey (Liu et al. in preparation) is expected to cover a larger 
solid angle especially at bright fluxes, and to be at least a fac- 
tor of two deeper. The exploitation of the Chandra and XMM 
archives is still far from being concluded, and we expect to see 
the number of high-z clusters (z > 1 ) to double in the next years, 
and in particular to see the increase of the sample of clusters at 
z > 1.5 which are being found only recently. However, a major 
breakthrough in the field of cluster surveys can be achieved only 
with a survey dedicated mission. 

The planned eROSITA satellite (Predehl et al. 2010; 
Cappelluti et al. 2010) and the proposed WFXT mission 
(Rosati et al. 2010; Murray et al. 2010) can provide a large num- 
ber of new detections and X-ray spectra for a large number of 
them. The eROSITA mission is expected fo fly in the near future, 
and it will finally provide an all-sky survey almost two orders of 
magnitude deeper than the last one performed by the ROSAT 
satellite (Voges et al. 1999), filling a gap of almost 20 years. 
Unfortunately, eROSITA will be confusion limited in the flux 
regime where the majority of the high-z clusters are expected. 
In addition, the low hard-band sensitivity (the effective area de- 
creases an order of magnitude 5 from 1 keV to 4 keV) hampers a 
measure of the temperature for hot clusters, limiting cosmologi- 
cal studies. 

The WFXT is the only X-ray mission optimized for efficient 
wide-area surveys, with a design which yields good angular res- 
olution and a constant image quality across a large (1 deg 2 ) field 
of view (see Rosati et al. 2010). These aspects, coupled with a 
large effective area, will provide not only a survey more than 
one order of magnitude deeper than eROSITA, but a direct mea- 
surement of global temperatures, density profiles and redshifts 
for a significant fraction of the cluster sample, thus allowing 
cosmological studies independently of optical and spectroscopic 
follow-up (Borgani et al. 2010). Simulations based on the orig- 
inal WFXT design (Tozzi et al. 2010) show that such a mission 
is able to match in depth, survey volume and angular resolu- 
tion surveys at other wavelengths which are planned for the next 
decade. 

9. Conclusions 

We present a new sample of X-ray selected groups and clusters 
of galaxies obtained with the X-ray Telescope (XRT) on board 
of the Swift satellite. We search for extended sources among 
336 GRB fields imaged with XRT with galactic latitude \b\ >20° 
(available in the archive as of April 2010). We identify extended 
sources with a simple criterion based on the measurement of the 

5 http://www.mpe.mpg.de/erosita/ 
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HPR within a box of 45 x 45 arcsec of the source image. We ap- 
ply a sharp threshold of 100 net soft counts within the extraction 
radius R ext to select reliable extended sources. Extensive simula- 
tions showed that our method, despite being very simple, provide 
us with an X-ray sample with an expected high completeness 
and a low contamination, also thanks to a careful a posteriori 
visual inspection. 

Our final group and cluster catalog consists of 72 X-ray 
sources. The sky coverage of the survey goes from the total 40 
deg 2 to 1 deg 2 at a flux limit of about 10~ 14 erg s _1 cm~ 2 . The 
corresponding logN-logS is in very good agreement with previ- 
ous deep surveys. 

We directly verified that there is no correlation between the 
position of the clusters and that of the GRBs. In other words, se- 
lecting our cluster sample from XRT GRB fields does not result 
in any spatial bias with respect to the GRB positions. 

A cross correlation with X-ray catalogs shows that only 9 
SXCS sources were previously identified in the X-ray band, 
none of them classified as extended. A search in optical 
databases (mostly based on SDSS data) allows us to find the 
counterparts of 20 clusters. In addition, 4 galaxies with redshift 
have been found to be within 7 arcsec from the X-ray centroid, 
and therefore they are considered as possible identification of 
the central galaxy of the group/cluster candidate. Overall, only 
20 sources are confirmed as clusters using data from the lit- 
erature, while 30 sources have some counterpart in the NASA 
Extragalactic Database, and, finally, 42 sources are newly iden- 
tified. 

We estimate that about one third of the sample is detected 
with a S/N high enough to allow the measure of the redshift from 
the X-ray spectral analysis (see Yu et al. 201 1), as we will show 
in a companion paper (Moretti et al. in preparation). Thanks to 
the quality of our X-ray data, and thanks to the synergy with 
other surveys, specifically the SDSS, we are able to provide a X- 
ray cluster sample with a well established completeness function 
down to a flux limit comparable to that of the deepest cluster sur- 
veys based on ROSAT data (RDCS, see Rosati et al. 1998), and 
with a comparable statistics. In addition, we expect to detect a 
few clusters with redshift z > 1 . Overall, the SXCS is expected to 
give a significant contribution in the field of X-ray clusters sur- 
veys also thanks to its peculiar properties of a low background 
and a constant PSF. These properties are two key requirements 
for future wide-area, X-ray surveys, as foreseen by proposed 
future missions aiming at bringing the X-ray sky to the same 
depth and richness of the optical and IR sky in the next decade. 
A deeper, extended release of the SXCS based on a new detec- 
tion algorithm tailored to XRT images is currently undergoing 
(Liu et al. in preparation). Catalog and data products of SXCS, 
constantly updated, are made avalilable to the public through the 
website http : //adlibitum. oats . inaf . it/sxcs. 
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SWJ124312+170451 SWJ093749+153540 SWJ092719+301342 SWJ090946+415713 




Fig. 16: Selection of SXCS cluster candidates imaged by the Sloan Digital Sky Survey (SDSS) with X-ray contours overlaid in 
green. The image size is 3x3 arcmin. The X-ray contours are computed on XRT images in the 0.5 - 2 keV energy band with the 
Chandra Ciao software, and the contours are drawn at respectively 9, 15, 30, 50, 90, 160 and 300 cr, where <x is given by the local 
background level. 
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Table 2: SXCS sources. 
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ra 
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Table 2: continued. 
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SWJ175640+332928 
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269. 16901 
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3.56 


32 
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12 
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SWJ181053+581527 


III. /2318 


58.25763 
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3.16 
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SWJ194004+782419 
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4.8± 0.7 


SWJ203723-440141 


309.34961 


-44.02815 
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SWJ215507+164725 
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SWJ215354+165348 
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57331 


5.66 
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SWJ222600-571248 


336.30238 
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J70X 


_ 1 
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33.6± 4.9 


SWJ222443-022031 


536.18286 


-2.3421 / 
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SWJ222953+194354 


337.47308 


19.73191 
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SWJ224207+233354 


340.53061 


23.56503 
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17 
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13. 1± 1.7 


SWJ230754-681505 


346.97769 


-68.25166 


82185 


2.75 


23 


147+ 


13 


10.6 


4.5+0.6 


SWJ230650-680401 


346.71140 


-68.06695 


82160 


2.75 


54 
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SWJ232248 +054809 


350.70193 


5.80263 


200154 


5.04 


39 


1549+ 


:42 


36.1 


20.8± 2.2 


SWJ232345-3 13048 


350.93924 


-31.51346 
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34 
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20 
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SWJ233518-662139 


353.82767 


-66.36105 


111461 
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30 


226± 17 


13.0 


5.1+ 0.7 


SWJ2336 17-3 13626 


354.07169 


-31.60723 


43325 


1.22 


54 


1400± 39 


35.1 


78.5+9.1 
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Notes. (1) SXCS Id; (2), (3) Coordinates of the X-ray centroid; (4) Exposure time corrected for vignetting; (5) Galactic absorption; (6) Extraction radius R ex , in arcsec; (7) Soft net counts wifhifij 
R exl ; (8) S/N; (9) Soft band flux corrected for Galactic absorption. j» 



Table 3: List of counterparts of SXCS sources. 



Name 


Cluster 


Distance (') 


redshift 


Galaxy 


Distance (') 


redshift 


X-ray 


Distance (') 


SWJ005500-385226 


EDCC 


0.952 




LCRS 


0.066 


0.164127 






SWJ02 1747-500322 














1AXG 


0.571 


SWJ023924-250504 








2DFGRS 


0.109 


0.1737 






SWJ035259-004342 


WHL 


0.102 


0.3251 (ph) 


SDSS 


0.09 


0.334925 






SWJ084749+133141 


WHL 


0.065 


0.363 (ph) 


SDSS 


0.061 


0.348628 






SWJ090946+415713 


WHL 


0.128 


0.14855 (ph) 












SWJ092729+301048 


WHL 


0.072 


0.365 (ph) 


SDSS 


0.073 


0.293 (ph) 


1AXG 


0.35 


SWJ092719+301342 


WHL 


0.223 


0.29975 (ph) 












SWJ093045+ 165931 








2MASX 


0.053 


0.177278 






SWJ093749+ 153540 


AMF 


0.0363 


0.2939 (ph) 












SWJ094816-131644 














CXO 


0.107 


SWJ101341+430655 








SDSS 


0.043 


0.448860 






SWJ105946+534809 


SDSS-C4-DR3 


0.065 


0.072 (ph) 


MCG 


0.066 


0.071136 






SWJ1 1581 1+452906 


AMF 


0.558 


0.4051 (ph) 












SWJ 123620+285905 


AMF 


0.272 


0.2305 (ph) 


2MASX 


0.053 


0.228573 






SWJ124312+170451 


NSCS 


0.73 


0.1424 (ph) 












SWJ 13 1300+080259 


WHL 


0.076 


0.5598 (ph) 












SWJ131522+164145 














1WGA 


0.558 


SWJ133055+420017 








2MASX 


0.068 


0.061154 






SWJ140637+274349 














1WGA 


1.083 


SWJ 140639+273546 


GMBCG 


0.211 


0.243 (ph) 












SWJ140728+274917 


Abell 


0.554 


0.17244 








1WGA 


0.796 


SWJ143646+275157 


AMF 


0.706 


0.2648 (ph) 












SWJ151324+305738 


CGCG 


0.092 


0.0717 












SWJ155742+353023 


MaxBCG 


0.547 


0.1549 


2MASX 


0.117 


0.158877 


RXJ/WGA 


0.095/0.21 


SWJ215507+164725 














CXO 


0.062 


SWJ222600-571248 


Abell 


0.268 


0.13 












SWJ2225 16-020827 














1WGA 


1.393 


SWJ232248+054809 


NSCS 


0.914 


0.45 (ph) 












SWJ233617-313626 


Abell 


0.63 


0.0623 (ph) 













Notes. (1) SXCS Id; (2) Cluster Optical catalog (3); distance from the X-ray centroid in arcmin; (4) cluster redshift; (5) Galaxy Optical catalog; (6) distance from the X-ray centroid in arcmin; (7) 
galaxy redshift; (8) X-ray catalog; (9) distance from the X-ray centroid in arcmin. 
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